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Higher-order pQCD corrections play an important role in the reproduction of data at 
high transverse momenta in the energy range 20 GeV < ^/s < 200 GeV. Recent calcu- 
lations of photon and pion production in pp collisions yield detailed information on the 
next-to-leading order contributions. However, the application of these results in proton- 
nucleus and nucleus-nucleus collisions is not straightforward. The study of nuclear effects 
requires a simplified understanding of the output of these computations. Here we summa- 
rize our analysis of recent calculations, aimed at handling the NLO results by introducing 
process and energy-dependent K factors. 

The main motivation of ultrarelativistic heavy-ion collision experiments is their role in testing quantum 
chromodynamics (QCD) and, in particular, their potential to display the predicted phase transition 
between hadron matter and the quark-gluon plasma (QGP) [1]. At higher transverse momenta {pt ^ 
3 GeV), QGP scenarios should be judged against the background of basic particle yields from primary 
nucleon-nucleon collisions, which can be described by perturbative QCD (pQCD) [2]. Thus, to verify the 
formation of the QGP state [3], we need precise knowledge about the 'partonic background' provided by 
nucleon-nucleon collisions. 

The pQCD-based description of particle production in proton-proton (pp) collisions is a challenging 
task. Using pQCD at relatively low energy further complicates the situation. Recent studies with 
higher order contributions [4-7] investigated the role of different QCD scales in detail. However, the 
application of these results as a pQCD reference is not an easy job for several reasons. First, a recent 
NLO pion calculation has difficulty in consistently reproducing the pion data [5]. Second, although still 
of relatively low order, these calculations are rather long and complicated. It would be beneficial to find 
a fast phenomenological shortcut with appropriate precision. Finally, a heavy ion collision is not a mere 
superposition of binary nucleon-nucleon collisions, but the incoming nucleons and the produced partons 
are subject to multiple scattering and to other medium effects (e.g. higher-twist [8], shadowing [9]) due to 
the presence of further nucleons or partons. The phenomenological description should have the flexibility 
of being able to include many-body effects. 

The first and widely-used step in creating such an effective pQCD description is the introduction of the 
so-called "ivT factor" , which, roughly speaking, renormalizes the leading term of the pQCD cross section 
and accounts for higher-order corrections in different processes [10,11]. The higher-order corrections are 
known to be large and process dependent [12-14]. Note that several different K factors exist in the 
literature. For example, from an experimental point of view, it is worthwhile to compare data to a given 
order of a pQCD calculation, and this is frequently expressed in terms of a if factor. Theoretically, the K 
factor can be thought of as the ratio of a full pQCD cross section to a cross section through given order, 
or as the ratio of calculations through different orders, or the K factor can be studied in the framework 
of a model. A study of the K factor in the framework of the Parton Cascade Model (PCM) was carried 
out recently [15]. 

Here we insist on defining the K factor within a given order of pQCD (specifically NLO) as the factor 
correcting the Born-approximation result in NLO to the full NLO result for the particular process. In 
any pQCD calculation the partonic cross section is convoluted with parton distribution functions, which 
describe some of the perturbatively incalculable soft physics. The NLO parton distribution functions are 
fitted to the data in the framework of an NLO calculation [16]. For consistency, it is important that 
parton distribution functions of the order of the cross-section calculation be used. The K factor can be 
introduced as: 

(Pp (Pp (Pp ' d^p 

In general, the K factor is energy and momentum dependent. However, in many cases it is approximated 
as a constant in a certain momentum region for a given energy. The K factor is traditionally thought to 
have values in the range of 1.5 - 2.5. 



In this paper wc summarize our recent results on the vahie and the behavior of the K factor in the 
energy region 20 GeV < -^i < 200 GeV. Within this interval, we concentrate on the lower energies, where 
proton-nucleus data are abundant. First we display our results on the "jet level" . We then investigate 
the modifications due to jet fragmentation and especially to pion production. We close with an overview 
of recent results on photon production. A forthcoming publication [17] will include comparisons to 
experimental data using K factors from the present paper. 

The computation of the partonic cross section through NLO is tedious. The calculations by the EUis- 
Kunszt-Soper (EKS) group [12,13] and by the Aversa-Chiappetta-Greco-Guillet group [18] are based on 
the matrix elements published in Ref. [19]. We use EKS's public FORTRAN-code [20] to calculate the 
Born and correction contributions in our jet-level analysis. For the parton distribution functions we 
use the MRST (central gluon) set [16], which incorporates a large amount of experimental information 
evaluated at the NLO level. 

The choice of the factorization/renormalization scale Q influences both the LO and the NLO jet cross 
sections. The typical value of this parameter varies around Pt/3 ^ Q ^ "ipx in standard calculations. 
Here we take Q = pt/2. (Naturally, pT refers to the transverse momentum of the jet in jet-level 
expressions, while it denotes the transverse momentum of the final hadron when observable hadron 
production is calculated.) It should be kept in mind that several complications arise in NLO which do 
not plague the LO calculation. Among these are the dependence of the NLO results on tlw^ jcrt-cone 
angle, R = \J (At;)^ + (AcE>)2 (with pseudorapidity r\ and azimuth <&) and the jet-separation parameter 
i?sep- Following EKS [20], the collinear singidarity, which also appears first in NLO, is treated here by 
choosing QcoiHnear = Pt- We are particularly interested in the dependence on R, and keep Rgep = 2i? 
in our calculations. We display two results. In Fig. 1 we use the 'optimal size' R = 0.7 for the jet cone, 
where the scale dependence on the factorization scale Q is minimal [21,22]. Fig. 2 displays results with 
a larger cone size, R = 1.0, where the bremsstrahlung contributions are enhanced, causing an increase of 
the K factor. In both figures we show K as a. function of pt for energies 20 GeV < ^/s < 200 GeV. We 
present calculated results for the px region where the accuracy of the code is better than 5%. The solid 
lines represent a fit to the calculated data points at the particular energy (linear fit in Fig.l, quadratic 
in Fig. 2). It can be seen in both figures that the magnitude and the transverse-momentum dependence 
of the K factor decreases with increasing energy. At ^/s = 200 GeV, the value varies less than 5% from 
Pt = 3 GeV to px = 10 GeV. In contrast, at yG — 20 GeV, the variation is on the order of 20% from 
Pt = 2.5 GeV to pr = 4 GeV. This is not unexpected, since higher-order terms should become less 
important with increasing energy. Comparing Fig.-s 1 and 2, we see that the K factor is larger if R 
is larger and the dependence of K on pt is more linear with the smaller value of R. This is due to 
the increased contribution of higher values of px with increasing cone size. For practical purposes, we 
parameterize the energy and momentum dependence of the K factor for i? = 0.7 as 
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For i?= 1.0 we find 
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These parameterizations are also shown on Fig.-s 1 and 2 (dashed lines). 

So far, we only discussed the K factor at the jet level. The above parameterization characterizes 
an average enhancement of jet production through NLO. To calculate observable hadron production via 
parton fragmentation in an NLO treatment, one may substitute the above Kjgf{s,pT) into the pQCD 
calculation leading to an approximate NLO-result: 
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where fa/p'^i^a- Q^) stands for the parton distribution function and D^j^^{zc,, Q'^) denotes the fragmen- 
tation function (both in NLO). The label c represents the fragmenting parton. The fragmentation scale is 
taken to be Q' = pt/2, while the factorization scale is kept at Q = pt,c/'^- Applying this approximation 



and using a set of fragmentation functions [23], one can easily determine e.g. pion and kaon production 
through NLO, and extract an appropriate K factor for (average) pions, K^{s,pt), defined as the ratio 
of the cross section in Eq. (4) to the cross section with Kj^t = 1- This is displayed in Fig. 3, together 
with a similarly- defined Kk{s,Pt), in the energy range 20 GeV < ^/s < 200 GeV. We obtained very 
similar pion and kaon K factors. One advantage of using Kjet is that a single approximate function can 
be applied in calculating the cross sections for several similarly-produced hadrons. 

The K factor at the hadron level is larger than its iet-lc;vel counterpart by 10 - 20%. This difference 
can be thought of as a 'backshift' in px, as some of the transverse momentum of the jet manifests itself 
in other particles after hadronization. Fig. 4. illustrates these results at the lowest energies considered, 
y/s = 24 GeV. Extracted Kjet (dotted lines) and (solid lines) as functions of pt are displayed for 
R = 0.7, 0.85, and 1. It can be seen that the pion K factor is approximately parallel to the corresponding 
jet-level K factor, shifted toward lower values of pT, as discussed above. 

Fig. 5 completes the picture. Here we present calculated K factor for photons, based on Ref. [4,24], 
using Eq. (1), in the same energy and transverse momentum range. It is important to note that this 
calculation includes the contribution from jet fragmentation [4]. We see that the Kj obtained in this 
manner shows a similar tendency to the K factor with R — 0.7 (Fig. 1). 

These results suggest that an averaged, approximate version of an NLO calculation represented by 
a K{s,pt) factor will be useful in heavy-ion collisions in this transverse-momentum range from SPS to 
RHIC energies. The above parameterizations (or similar ones at different R) can be utilized together with 
the Born calculation in this region. As an application, and as an improvement on our earlier work [25] , 
we intend to use these K factors in the energy range 20 GeV < ^/s < 60 GeV, at transverse momenta 
3 GeV <Pt<7 GeV [17]. At RHIC energies (V^ = 140 - 200 GeV) we find an approximately constant 
K factor. The numerical value with Q = pt/2 and R = 0.7 is K 1.2. 

In conclusion, we have developed a simplified phenomenological description of recent NLO results for 
particle production in pp collisions in terms of energy and transverse- momentum dependent K factors. 
Since the complexity of proton-nucleus and nucleus- nucleus collisions presents its own challenges, we 
anticipate that our parameterization of the K factor will find useful applications in the analysis of 
forthcoming data. 
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Fig. 1. Energy and transverse-momentum dependent K factor, Kjet{s,PT) with Q = pt/^, 
R — 0.7, Rsep — 2-R. Dashed Hnes represent Eq. (2). 
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Fig. 2. ii' factor, Kjet{s,PT) with i? = 1.0 (Q = Pt/2 and i?sep = 2R). Dashed hnes correspond 
to Eq. (3). 
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Fig. 3. K factor for pfons, Kj^^SjPt) (solid line) and for kaons, Kk{s,pt) (dashed), after 
hadronization of jets calculated with R = 0.7 {Q = pt/2, Rsep = 2i?). 
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Fig. 5. K factor for photons, Kj{s,pt) at energies y/s = 20 — 200 GeV. 



